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Abstract 

The purpose of this report is presentation of the main modifica- 
tions of the standard Kibble-Zurek formalism caused by the existence 
of unperfections in the system. We know that the distribution of 
kinks created during a second order phase transition in pure systems 
is determined solely by the correlation length at freeze-out time. The 
correlation length at that instant of time intuitively describes the size 
of the defect and therefore the number density of defects is limited 
by the possibility of holding kinks in a unit volume. On the other 
hand if the system is populated by the impurities then kinks emerge 
mainly in knots of the force distribution which correspond to extremes 
of the impurity potential i.e. positions of imperfections. The purpose 
of this report is to show that, due to existence of the strong gradients 
of the impurity potential, kinks can be created mainly in the close 
vicinity of the impurities. It seems that this simple mechanism can 
be responsible for occurrence, in the number density formula, the ad- 
ditional length scale describing the impurity distribution. We know 
that in pure systems, as a consequence of kink-antikink annihilation, 
the number density of kinks decrease in time. In contradiction to pure 
systems, kinks produced in the systems populated by impurities could 
be confined by the impurity centers and therefore they may not dis- 
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appear from the system and may remain above the level established 
by thermal nucleation of pairs. 

1 Introduction 

Last years topological defects attract attention of many researchers. The 
motivation of these studies comes from the fact that they can be seen as 
macroscopic manifestations of underlying physical processes. On the other 
hand they can help to study the nature of critical dynamics. 

The theory describing the dynamics of the second order phase transition 
was proposed by Kibble and Zurek [1]. The key point of the Kibble-Zurek 
mechanism is an observation that the order parameter evolves adiabatically 
through a sequence of nearly equilibrium configurations up to the freeze-in 
time. At that instant the system loses capacity to respond for the changes 
of the external parameters. From that time to the freeze-out time the field 
configuration remains almost unchanged. The dynamical evolution restarts 
below the critical temperature at freeze-out time. At that instant the sys- 
tem regains ability to respond for the changes of the external parameters 
but it is too late to undo non-trivial arrangements of the order parameter 
from above the critical point. This paradigm works for the overdamped and 
underdamped systems as well. The main prediction of this scenario is the 
dependence of the number density of produced defects on correlation length 
n ~ £~ d at freeze-out time or its dependence on quench time n ~ r~ d ^ 4 , 
where d denotes the number of space dimensions. This scenario was well 
verified in a series of numerical experiments [2]. 

The defect network density obtained at freeze-out time is an initial condi- 
tion for dynamics which is determined by the defect - antidefect interactions. 
Due to annihilation of defects and antidefects the initial density of defect 
network first is quickly reduced in time and then is stabilized on the level de- 
termined by the Boltzman factor which describes the probability of thermal 
nucleation of the kink-antikink pairs [3]. 

In real life experiments researchers use mainly stable coherent quantum 
systems. So far, experiments were performed in Helium-3 on symmetric phase 
3 He-B which is more simple to experimental and theoretical treatment. The 
results of experimental studies confirms the dependence of the number den- 
sity of produced vortices on quench time [4]. Importance of studies of the 
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transitions in Helium-3, follows from the fact that due to nontrivial struc- 
ture of the order parameter it allows for experimental verifications of ideas 
concerning the structure of vacuum of the quantum field theory. 

Researchers performed also experiments on liquid crystals. The creation 
of disclinations of different types produced during a quench from disordered 
to nematic phase in liquid crystals was examined and the results were to 
some degree consistent with the Kibble-Zurek predictions [5]. 

There were also more controversial experiments made in superfluid helium- 
4 where almost no vortices of topological origin were observed [6]. 

Lately there are also attempts to study the creation of vortices in optically 
cooled alkali atom clouds during formation of the Bose-Einsten condensate. 

In this report I would like to concentrate on influence of impurities on 
creation of topological defects. It is difficult to imagine free of unperfections 
liquid crystal or even superconductor. The population of the superconduc- 
tors and liquid crystals by the impurities and admixtures seems to be an 
inevitable outcome of their preparation. On the other hand quantum liquids 
are one of the purest substances in the nature. Although the solubility of 
foreign materials in liquid helium is almost zero there exists some artificial 
techniques, like aerogel technique [7], which allow to introduce impurities 
even into quantum liquid. 

Prevailed part of the obtained hitherto results concern homogenous medium. 
On the other hand, the presence of the impurities can significantly change 
properties of the system. 

This report aims in presentation of the main modifications of the stan- 
dard Kibble-Zurek formalism caused by the existence of unperfections in the 
system. 

2 An influence of inhomogenities on produc- 
tion of kinks 

First let us recall K-Z formalism applied to description of homogenous sys- 
tems. A pure, overdamped 4 system is described by the following equation 
of motion 

7<9 t 0(t, x) = dl<f)(t, x) - a{t)(f)(t, x) - A0 3 (t, x) + rj(t, x), (1) 
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where rj(t, x) is a temperature white gaussian noise defined by the correlators 



{V(t,x)) = 0, 



(rj(t : x)rj(t' : x')) 
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5(x-x')5(t-t'). 



(2) 



The explicit time dependence of the "mass" parameter a(t) allows for mod- 
elling the phase transition in the system. Depending on the sign of this 
parameter we can find a system in the phase with one trivial or two nontriv- 
ial ground states. 

The number of kinks produced during a phase transition at freeze-out 
time is calculated from the Liu-Mazenko-Halperin formula [8] 



In fact this formula can be expressed with the use of the power spectrum 
which is defined by the equal time correlator of the order parameter. The 
cut-off in this formula separates the stable from unstable modes of the system. 
We integrate only over unstable modes because only they can grow to form 
stable kink structures. As a result of this calculation one could obtain a 
Kibble-Zurek critical exponent 1/4 which describes the dependence of the 
number density of produced kinks on quench time r 



In this context usually a linear quench is presumed. The critical exponent 
depends only on number of spatial dimensions. 

On the other hand, many physical systems are dense populated by imper- 
fections of different types. The presence of impurities in the system can be 
taken into account by introducing a deterministic force distribution V(t,x) 
into the equation of motion 



7<9 t 0(t, x) = d 2 x (p(t, x) - a(t)(f)(t, x) - A0 3 (t, x) + r)(t, x) + V(t, x). (5) 



First we reconsider the way of counting zeros of the order parameter. The 
number density of produced zeros can be defined as a ratio of zeros located 
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(3) 
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in some interval of space to the length of this interval. On the other hand 
the number of zeros can be calculated as a sum of arbitrary quantity divided 
by itself over all points where the scalar field disappears. One of the possible 
choices of this quantity is 0' 

We know that integration of the delta function can be replaced by the sum- 
mation over zeros of delta argument 



dxf(x)5[g(x)] = ^2 | ^ ^| where Xi are defined by the equation g(xi) = 0. 



If we identify in this lemma function / with 0' and g with <p then we replace 
the sum over all zeros by the average of some integral 
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[X) = l™2L { Ll <W(M)IW(M)]>- (7) 



In zero L limit this expression reduces to the average over all realizations of 
the noise of some combination of the delta and sign functions 

n{x) = (sign[(f)'(t, x)](j>'(t, x)6[(f>(t, x)}). (8) 

In the next step we use integral representations of the signum and delta 
functions and then divide the scalar field on two components <p(t, x) = 
ip(t,x) +u(t,x). One which carries deterministic part of the evolution u 
and the second ip which carries the stochastic part of the evolution of the 
scalar field 0. We known that above critical point the field fluctuates around 
trivial ground state. If the amplitude of the noise is small i.e. if the temper- 
ature of the system is low then fluctuations of the field are also small and 
its value is close to zero. As the field fluctuate around zero value its aver- 
age magnitude is small and therefore cubic term in the equation of motion 
is negligible. As the system evolves through a sequence of almost equilib- 
rium states this remains true up to freeze-in time. Identification of those two 
components in linear approximation is straightforward. 

Then we use theorems which allow to replace n-th order correlators by the 
average and correlator of second order and therefore we are able to replace 
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the averages of some functions of random variable by the functions of the 
second order correlator of this variable. 

Finally we obtain a formula which in the absence of impurities reduces 
to the well known Halperin-Liu-Mazenko formula (3) 



, \ 1 
n{t, x) = — A 



(ib 12 ) " 2 . " ,2 . u' " 2 . / u ' 

xr '. e W 1 ) ^W 1 ) + —— e *iF>Erf[ 1, (9) 



M (^ 2 ) v / 2tt(V' 2 ) V\/ 2 ^ /2 ) 

where Erf is the error function. The analytical and numerical studies of 
the kink distribution shows that kinks are created mainly in the vicinity of 
knots of the force distribution which corresponds to extremes of the impurity 
potential [9]. 

In most of the physical systems the shape and the distribution of the 
impurities is random and therefore we have to allow the force distribution 
to be a random type with some length scale which characterizes the average 
distance between impurities. This time the equation of motion contains two 
random forces. First represents thermal fluctuations in the system 77 and the 
second one which describes distribution of impurities V. 

The angle bracket represents the average with respect to all realizations 
of the thermal noise and the new bracket {...} represents an average with 
respect to all possible distributions of the impurities in the system. 

Equation of motion in this new situation, at first sight, seems to be iden- 
tical with equation (5) 

7<9 t 0(t, x) = d 2 x <f){t, x) - a(t)(f)(t, x) - \<f) 3 (t, x) + r)(t, x) + V(t, x). (10) 

In fact 7 in this equation is not a simple constant but it is an integral oper- 
ator. The existence of this term is an inevitable if we restrict our studies to 
stationary processes. The explicit form of this operator is the following: 

rt /"CO 

7 <9 t 0(t,a;)= / dt' d 3 x'-f(t, t'; x - x )d t 4(t' , x ). 

J to J — OO 

If the impurity force distribution T> has a form of the gaussian white noise 
then ry(t,t';x — x') = 7#(t — t')S(x — x'), and this integral reduces to the 
damping constant multiplied by the time derivative of the order parameter. 
In generic situation we expect a dependence of V correlators on some length 
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scale (e.g. average distance between impurity centers) and therefore the force 
distribution is defined as follows 



{V(t,x)} = 0, 
{V(t,x)V(t',x')} = \w{\x-x'\)5{t-t'). 



(11) 



To better understand the source of this complication let us consider a 
simple mechanical analogy which is a Brownian motion theory. 

The erratic motion of a Brownian particle is caused by collisions with the 
molecules of the fluid in which it moves. These collisions allow an exchange 
of the energy between the fluid and the Brownian particle. If the Brownian 
particle is much more massive than the molecules of the fluid then the in- 
fluence of the molecules on the observed particle can be approximated by a 
Gaussian white noise rjcit): 



where x(t) is the position of the Brownian particle. The generalization of 
Brownian motion theory to the random motion of a particle which is not 
necessarily heavier than the molecules of the fluid was proposed by Kubo 
[10]. In this case the time scale of molecular motion is no longer very much 
shorter than that of the motion of the particle under observation, so that 
the random force r)(t) can not be of Gaussian type. To describe an influence 
of the molecules on the observed particle we have to introduce a color noise 
characterized by some time scale. This time scale may describe the average 
time interval between two subsequent collisions of the molecules with the 
observed particle. In addition, if we consider a stationary process we have to 
abandon the assumption of a constant friction and to introduce generally a 
frequency-dependent friction 



In case of 4 model the distribution of impurities is not generally described by 
the white gaussian noise. In generic situation the distribution of impurities 
is characterized by some length scale which describes the average separation 
o impurity centers and therefore, in similar way as it was in case of Brownian 
particle, we introduce retardation to </> 4 model. 



mx(t) + jx(t) = i] G (t), 



(12) 




(13) 
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Next we have to make further generalization of the Liu-Mazenko-Halperin 
formula. This generalization is achieved by averaging the formula (9) with 
respect to possible distributions of the impurity centers 



n = {n(t, x)} 



(V>' 2 ) + {W 2 } 



(14) 



The further understanding can be made for particular choice of the noise 
amplitude. The most representative one is Ornstein-Uhlenbeck amplitude 

\x\ 

W(\x\) = Ae~~ which interpolates between constant distribution and gaus- 
sian white noise. In this model the number density of produced kinks depends 
on quench time and on characteristic length scale of the impurity distribution 
as well [11] 



n 



1 



0.434- + 0.34 



L 2 



(15) 



where b 



-r and c 



1.816 + 0.83c ' 
-^=. Let us notice that in case of week imperfections 
the usual scaling is recovered 

n~\. (16) 

T4 



3 Remarks 

We show that there are two components which determine the number density 
of kinks produced in the inhomogenous system during the second order phase 
transition. First component follows from the K-Z formalism for pure systems 
and is determined by the quench time. The second component is determined 
by the characteristic length which describes the distribution of impurities in 
the system. Due to existence of impurities and admixtures the kinks are 
created mainly in the knots of the impurity force distribution. One could 
even find exact solutions which describe the kinks confined by some partic- 
ular impurity potentials [12]. An examples of those solutions are squeezed 
kink or squeezed anti-kink. One could also check linear stability of those 
solutions. The other solution obtained for some particular force distribution 
is a static kink-antikink solution which illustrates that, in contradiction to 
pure systems, the configuration of this type have no tendency to annihilate. 
This is also a reason why at late time after transition in the inhomogenous 
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system we still can find the number density of kinks substantially larger than 
estimated from the probability of thermal nucleation of kink-antikink pairs. 
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